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Abstract 

We investigate metastability in the two dimensional Ising model in a square 
with free boundary conditions at low temperatures. Starting with all spins down 
in a small positive magnetic field, we show that the exit from this metastable 
phase occurs via the nucleation of a critical droplet in one of the four corners of 
the system. We compute the lifetime of the metastable phase analytically in the 
limit T — ^ 0, /i ^ and via Monte Carlo simulations at fixed values of T and h 
and find good agreement. This system models the effects of boundary domains in 
magnetic storage systems exiting from a metastable phase when a small external 
field is applied. 
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1. Introduction. 

Metastability is observed in many different systems close to a first order phase tran- 
sition. It is a dynamical phenomenon (see, for instance, [PL]) not included in the Gibbsian 
formalism, which is so successful for the description of stable equilibrium states [LR,I]. 
The development of a full theory of metastability is desirable for its intrinsic as well as 
for its experimental and technological interest and it also poses challenging mathematical 
problems (see [CG0V,0S1,0S2]). 

The metastable behavior of the nearest neighbor two dimensional Ising model for 
large finite volumes and small magnetic fields was analyzed in [NS1,NS2] in the zero tem- 
perature limit in the framework of the "pathwise approach" introduced in [CGOV]. In 
[SI] R. Schonmann, using arguments based on reversibility, described in detail the typical 
escape paths from the metastable to the stable regime. Other regimes, very interesting 
from the physical point of view and mathematically much more complicated (finite tem- 
perature, infinite lattice and zero magnetic field), are considered in [S2] and [SS]. The finite 
temperature case has also been widely studied by Monte Carlo methods, for instance in 
[B,BM,BS,TM1]; a complete and clear description of these numerical results can be found 
in [RTMS]. This case has also been studied by means of transfer- matrix and constrained- 
transfer-matrix methods in [PS1,PS2,GRN]. 

In the same asymptotic regime as in [NSl], different Ising-like hamiltonians have been 
considered in [K01,K02,N0] and the three dimensional nearest neighbor Ising model has 
been studied in [EC] . In [CO] the Blume-Capel model has been studied (see also [FGRN] 
for a study of a version of this model with weak long range interactions) and in [0S1,0S2] 
the problem of metastability has been investigated in a more general case. 

All the above works have been carried out for systems with periodic boundary con- 
ditions; in [RKLRN] Ising model has been studied in the case of semiperiodic boundary 
conditions. In this note we study the case of a finite lattice with free boundary conditions 
at low temperatures and small magnetic fields using both rigorous analysis and Monte 
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Carlo simulations. 

The case of free boundary conditions is of some technological interest: during the 
recording process on magnetic tapes different parts of the magnetized material consisting 
of fine magnetic particles are exposed to different magnetic fields, resulting in domains 
with different orientation of the magnetization. In order to be used as storage devices, 
these materials must be able to retain their magnetization for long periods in weak ar- 
bitrarily oriented magnetic fields. The study of the escape from a metastable phase in a 
periodic system neglects the effects of boundary domains. Such effects are modeled here 
by considering a lattice with free boundary conditions. 

We find that although the main features of the nucleation of the stable phase are not 
changed, some interesting new aspects arise: in particular one can say a priori where the 
nucleation of the stable phase will start, that is where the critical droplet will show up. 

The model and results are described in Section 2 and analyzed in Section 3 and 4. 
Section 5 is devoted to some brief conclusions. 

2. The model and the results. 

Let us consider a two dimensional Ising model defined on a finite square A = 
{!,..., M}^ C 1? with free boundary conditions. The space of configurations is denoted 
by f2 = {—1, +1}^ and to each configuration cr e is associated the energy 

^(^) = E ^(^)^(^) - ^ E ^(^) ^(^) = ±1 ' (2-1) 

<x,y> x€A 

where the first sum runs over all the pairs of nearest neighbors in A and J,h > 0. The 
equilibrium states are described by the Gibbs measure 

where /? is the inverse temperature. 
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The time evolution of the model is given by the Metropolis algorithm: given a 
configuration cr at time t, we pick a site x at random and then change cr{x) to —cr{x) with 
probability 1 if AH is < and exp{-(3 AH) if AH > 0. 

It is easy to show that this dynamics is reversible with respect to the measure (2.2), 
hence, the unique invariant measure of the process is the equilibrium Gibbs measure. The 
problem we want to study is the way in which a system approaches the equilibrium state 
when it is prepared in the configuration with all the spins equal to minus one (ctq = —1) 
and the magnetic field h is chosen positive but small with respect to the coupling constant 
J {j < 1), while P is very large. 

In [NS1,RTMS,S1,TM1] this problem was studied for periodic boundary conditions: 
it was shown that for /? sufficiently large and h small enough, depending on the size of the 
lattice, the system shows metastable behavior. This means that the system spends a long 
time Tp^p ~ exp(/3^^) in a phase with negative magnetization performing random wan- 
derings near the configuration —1. These wanderings are characterized by the formation 
of small droplets of pluses inside the sea of minuses which disappear quickly; their typical 
life time r^^ is much shorter than the lifetime of the metastable state (r^^ <^ Tv,I3)- 

A droplet of pluses will tend to grow however if it is large enough, i. e. when its 
dimension is larger than a critical length, I* = [^] + 1, where [a] is the integer part of the 
real number a. The exit from the metastable phase is achieved, then, when a sufficiently 
large droplet shows up somewhere in the lattice: this droplet is called protocritical and 
in the limit /? — > oo it is a square droplet with sides I*. When this protocritical droplet 
appears, it grows and covers the whole lattice in a time ^ which is very small compared 
to the lifetime of the metastable phase. It has been shown that /j r^^ -C 

In our case of free boundary conditions we show, rigorously in the limit /? — > 00 and 
via Monte Carlo simulations for /3 large, that the system exhibits metastability and that 
the lifetime of this metastable phase is rg ~ exp(/3^). This exit time is much smaller 
than in the case of periodic boundary conditions. This is a consequence of the fact that 
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the tendency of a droplet to grow is favoured when such a droplet has one of its sides on 
the boundary of the domain or at a distance one from it. Indeed, we have to introduce 
two different critical lengths, Ai and A2, which refer respectively to droplets close to and 
far from the boundary. In the limit /3 — > 00 we find 

J 2J 
Ai = [;^] + 1 and A2 = [-^] + 1 , (2.3) 

with J and h fixed such that ^ > 1- 

Finally, we remark that the exit from the metastable phase occurs through a critical 
seed which appears in one of the four corners of the lattice. It grows in a time ~ 
exp{P{J — h)) to cover the whole domain. Hence, in this case the position of the nucleation 
seed in the lattice can be predicted a priori. In [RKLRN] it was observed that in the case 
of semiperiodic boundary conditions the critical droplet can show up on one of the two 
sides where the periodic boundary conditions are not imposed. 



3. Numerical results. 

In this section we describe some simulation results obtained using the Metropolis 
algorithm for an M by M square at low temperature and small magnetic fields; the typical 
values which have been used in the numerical experiments are P > 2, h < 0.5, and J = 1. 
This range of parameters is different from those considered, for instance, in [RTMS,TM1] 
and references therein, e.g. in [RTMS] they considered the case J = 1, /? = 1.102 and h 
varying approximately in the range 0.04 < h < 0.9. 

To obtain a numerical estimate of the critical lengths Ai, A2 we fixed (3 and h and 
prepared the system in the configuration a = —1 except for a plus droplet of size I in 
the corner of the square; the size of the square was chosen M = I + 5. We considered 
decreasing values of I and in each experiment we found the smallest value of I such that 
the droplet grew. We took this as an estimate of the critical length Ai. The length A2 was 
measured in a similar way; the droplet now being placed at the center of the lattice. 
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Fig. 1 presents an average over 60 independent determinations of Ai and A2 with 
inverse temperature /? = 10, and /? = 6 respectively. The sohd and dashed hnes represent 
the theoretical values (2.3), valid in the limit /3 ^ 00, while the black circles and black 
squares are, respectively, the numerical estimates of the critical lengths Ai and A2. The 
error bars have been omitted because the statistical errors, evaluated as the empirical 
standard deviation over the square root of the number of experiments, were found to be 
very small. The agreement between the numerical result and the theoretical prediction is 
very good. 

We remark that the number of Monte Carlo steps per site (MCS) one has to wait in 
order to see the growth of the droplet greatly increases as the value of the magnetic field 
h is decreased; indeed this time, for /? large, is approximately given by exp/5(J — h) and 
exp P{2J — h) respectively in the case of a droplet close to or far from the boudary of the 
domain. The smallest magnetic field we have considered is = 0.02: in this case we have 
set M = 56 and the number of MCS needed to see the shrinking or growth of the droplet 
was approximatively 10^. 

In Fig. 2 we have plotted the magnetization per spin of the whole box tuq and the 
magnetization per spin mi, m2 mz and 777.4 evaluated in four square boxes of side Ai placed 
at the four corners of the lattice, as functions of the number of iterations (time) in a single 
history of the system obtained after preparing the system in the starting configuration — 1. 
In both pictures the solid line represents the magnetization of the whole lattice mo, while 
the other lines refer to mi, m2 ma and rriA. The top graph in Fig. 2 refers to the case 
where P = 3, h = 0.24 (Ai = 5) and M = 16; the one below has been obtained with /? = 2, 
h = 0.14 (Ai = 8) and M = 32. 

In both cases it is clear that the system stays for a long time (about 10^ MCS) in the 
configuration —1: the thermal fluctuations, visible for /3 = 2, are negligible for /3 = 3. After 
this long time the magnetization in one of the corner boxes flips to one (nuclcation of the 
protocritical droplet); once this rare event has happened, all the other magnetizations start 
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to grow and quickly reach the value +1; that is the system quickly reaches the equilibrium 
state. 

Finally, we have evaluated the lifetime of the metastable state. In [RTMS,MT] and 
references therein the lifetime of the metastable state has been estimated at temperature 
T = 0.8Tc, where Tc is the Onsager critical temperature, corresponding to our (3 = 1.102. 
As the magnetic field is varied, four different regimes are detected (see [RTMS] for a 
detailed description of these different regimes); at low values of the magnetic field they 
find that the logarithm of the lifetime of the metastable state is a linear function of the 
inverse of the magnetic field. They call this the "single-droplet region" , meaning that the 
nucleation of the stable state is achieved via the formation of a single cricical droplet, see 
also [RKLRN,TM2] 

In Fig. 3 we have plotted our numerical measurements of the escape time versus 
the inverse of the magnetic field. It is clear that when the temperature is decreased the 
numerical results approach the theoretical value 

ilogr.^i (3.1) 

which by Theorem 1.1 is valid in the limit of zero temperature. The numerical data have 
been fitted with a linear function ^ logrg = to/j^ + n^j; the values of m/j and n^j are listed 
in Table 1. It is clear that the trend with increasing /3 is correct, although the value (3.1) 
would be reached only for much larger values of (3. Such simulation would require very 
long computer time. 

The results in Fig. 3 have been obtained in the case M = 32 and are the average of 60 
different histories. We have checked that very similar results are obtained if one considers 
larger domains (for instance M = 64, 128) but have not performed extensive statistics in 
these situations because the behavior of at small values of the magnetic field and very 
low temperatures does not depend on the size of the lattice (this is confirmed by the results 
in [RTMS], see Fig. 2 there). Finally, we remark that the linearity of the logarithm of Tfs 
is lost when 1/h is small enough, this is because for h large enough the system is in the 
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multi-droplet regime. 

4. Rigorous results. 

We list some definitions and notations which are necessary to discuss our rigorous 
results: 

1. The spin configurations a, ry G are called nearest neighbor configurations iff 3a; G A 
such that rj = a^, where is the configuration obtained by fiipping in a the spin at site 
X. A path is a sequence of configurations aoai...an such that for i = 1, ...,n — 1, crj_i 
and CTj are nearest neighbors. A path aoai...an is called downhill iS H{ai+i) < H{ai) for 
i — 0,1, ...n — 1. 

2. Given A <Z Q, rj & Q we define the hitting time 

rl = mf{t > : a'^ e A} 

3. A local minimum of the Hamiltonian is a configuration a such that one has 

iy(c7^) > H{a) Vx G A . (4.2) 

A local minimum will also be called a stable configuration because starting from it the 
system will not move for a time which is exponentially long in the inverse temperature /3. 

4. The set of all the local minima is denoted by C ri. 

5. Given G Al, we consider the process starting from r] and say 

rj subcritical <^=^ lim P{t^i < t1^) = 1 

/3— >oo - - 

r] supercritical <^=^ lim P{t1i < t^i) = 1 

/3— »oo - - 

where P{X) is the probability of the event X. 

6. Given u G we define its basin of attraction 

B{a) = {rj & Q : all downhill paths starting from rj end in a} . (4-4) 
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We observe that a downhill path oj = aoai...an necessarily ends in a local minimum if 

^ _ max^gn H{r]) - min^gn H{r]) 
h 

7. Given ^ C il, ^ is connected iff Vr/, a & Q 3 a path uj <Z Q starting from a and ending in 
77; we will say that this path connects a to 77. 

8. Given a connected set ^ C il, we will call the boundary of Q the set 

dg = {r]&Vt: r]^g,3xeA: r/* G ^} . (4.5) 

9. Given two rectangles i?i and R2 on the dual lattice A + (|, |), we say that i?i and R2 
are interacting rectangles iff Ri and R2 intersect or are separated by one lattice spacing. 
If two such rectangles have only two corners at distance one, then they are considered not 
interacting. 

Let h > and start the system from the configuration ctq = — 1 . We will now describe 
how the stable configuration +1 is approached and we will evaluate how long the system 
remains in the metastable phase. To do this we first characterize the local minima of the 
hamiltonian (2.1). 

Lemma 1.1 

Let us consider model (2.1) with J > h> and M > 2. Then a G iff a(x) = -1 Vx G A 
except for sites which are inside some rectangles laying on the dual lattice 

A + (|, |) such that Vz, j = 1, ...,n and i ^ j 
i) Ri and Rj are not interacting 
ii) Ri has sides longer than two 

Hi) Ri cannot have one of its sides on the "border" of A and one of its two other sides 
perpendicular to this one at a distance one from the border (see Fig. 4). 

Remarks. 

1. We denote by 7^(/i, h) with 2 < /i, /2 < M the set of all configurations with all spins 
— 1 except for those inside a rectangle with sides li and I2 and such that they are local 
minima. 
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2. Given a local minimum cr e TZ{li, I2) we denote it by Ri^m where 

l = mm{li,l2} and m = max{ii,Z2} • (4-6) 

3. We denote by 7?. C the set of all local minima containing only one rectangle of pluses. 

Now, we state under which conditions a local minimum is subcritical, that is we 
describe the evolution of the system starting from a stable configuration. Consider a local 
minimum Ri^rn G ^• 

Lemma 1.2 

When each side of Ri,m is at least at distance two from the border of the lattice, then 
given e > one has 

i) I < X2 =^ Ri,m is subcritical and 

ii) I > A2 =^ Ri,m is supercritical and 

When at least one of the sides of Ri^m is at distance one from the border of the square or 
is laying on it, then given e > one has 

i) I < Xi, m < M — 1 =^ Ri,m is subcritical and 

p(^^f3il-l)h-f3e ^ ^Ri^.m ^ g/3(i-l)/i+/3£-) ^4_g^ 

ii) I > Ai =^ Ri,m is supercritical and 
Remarks. 
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1. In the second case, if m = M then the local minimum is supercritical, no matter how 
long the smallest side I is. 

2. If one considers a local minimum with more than one rectangle, then this local minimum 
is subcritical iff all its rectangles are "subcritical" (abuse of language). 

3. It is possible to prove a stronger version of Lemma 1.2 that gives rise to a more detailed 
description of the contraction or the growth of the droplet (see Lemma 3 and 4 in [KOI], 
Theorem 1 in [NSl]). 

The proof of this lemma is the standard one, see, for instance, [NS1,K01,K02] or 
the proof of Proposition 4.1 in [CO]. One has to consider the basin of attraction of the 
local minimum Ri^m and work out the minimum of the energy on its boundary. Once this 
has been done, everything follows via the general arguments in Proposition 3.7 in [OSl]. 
The difference with respect to the case of periodic boundary conditions is that one must 
take into account situations like the one depicted in Fig. 5, that are absent in the Ising 
model with periodic boundary conditions. 

We can formulate, now, the theorem which describes the exit from the metastable 
phase: this theorem states that with high probability the system will visit a particular 
configuration V before reaching +1 and that the exit time is dominated by the time the 
system needs to reach V. 

This protocritical configuration V is such that all the spin are minuses except those 
in the union of a Ai x (Ai — 1) rectangle, one of whose corners coincides with one of the 
corner of the domain, and a unit square laying on the border of A and touching the longer 
side of the rectangle (see Fig. 6). 

Setting 

F = H{V) - H{-r) = J + 2JAi - h{Xl - Ai + 1) ; (4.11) 

we consider the process at starting from —1 and define t_i as the last time the configuration 
was —1 before it became +1, 

f_i = sup{t < T+i : at = -1}. (4.12) 
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Similarly we call r the first time after t_i when at = V, 



T-p — mi{t > T-i : (Tt — V} . 



(4.13) 



Finally, we can state the following theorem: 



Theorem 1.1 




1 



From Lemma 1.2 and Theorem 1.1 we have a rather accurate description of the 
system in the metastable phase: starting from —1 the system will spend a lot of time 
"close" to this configuration; sometimes small droplet of pluses appear, but the system 
quickly goes back to —1. Only after a long time, compared to the time these fluctuations 
need, the system will nucleate the protocritical droplet V and it will then reach the stable 
phase in a relatively short time. 

It is possible to give a more detailed description of the first excursion from —1 to +1; 
one could state a result similar to the one in [SI] (see also Theorem 3 in [KOI]). We do not 
enter in the details of this construction for the case of free boundary conditions, because 
no relevant difference appears with respect to the case of periodic boundary conditions. 
Roughly speaking one can say that that the system, during the excursion from —1 to +1, 
will follow a rather well-defined sequence of configurations made up of growing rectangular, 
almost square, droplets located at one of the four corners of A. 

The proof of Theorem 1.1 is now sketched: as in [KOI] we can define a set A satisfying 
some properties which will be listed below. The construction of the set A is exactly as in 
[KOI], except for the fact that whenever there is a rectangle with one side on the "border" 
and another side at distance one from the border of A, the rectangle is enlarged so that 
the latter side also touches the border (see e.g. Fig. 4). 

The relevant properties of the set A are the following: 
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i) A is connected; — 1 e ^ and +1 A. 

ii) There exists a path oj connecting —1 with V contained in A such that 

H{a) < H{V) V(7 e a;, (7 7^ P . 
There exists a path u' connecting V with +1 contained in A'^ such that 

H{a) < H{r) Na^uj\ ai^V . 
in) The minimal energy in dA is attained for the "protocritical" configuration; namely 

min(iy(c7) - Hi-Vj) = HiV) - H(-l) = T (4.14) 



and 



mm {H{a) - H{V)) > . (4.15) 

(TeA\{V} 



iv) With probability greater than zero, uniformly in P, the system starting from V will 
reach +1 before visiting —1; namely, given e > 

P(rri < r^) > e-^^ 

and 

Using properties i) — iv) and Propositions 3.4, 3.7 in [OSl] we get Theorem 1.1. 

5. Concluding Remarks. 

In this paper we have studied metastability in the two-dimensional Ising model on 
an M by M square with free boundary conditions rigorously in the limit — > oo and 
via Monte Carlo simulations at finite temperatures. We found good agreement between 
the theoretical predictions and the simulations and for a large range of h and low enough 
temperatures. The qualitative agreement persists even above one half of Tg. 
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Compared to periodic boundary conditions there are two relevant differences: i) the 
critical length of the droplet and hence the life time of the metastable phase is much 
shorter; ii) the protocritical droplet is always at one of the four corners of the square. 

It is clear that our analysis applies equally well to a rectangular domain with suf- 
ficiently long sides. In fact the basic approach carries over, in principle, to a general 
domain with general boundary conditions. The protocritical domain will always form, 
when /3 — > cxD at the place (or places) where the energy cost, H{V) — H{a-i), is minimal. 
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Table 1 









1.00 


0.216 


1.816 


1.25 


0.312 


1.173 


1.50 


0.483 


0.669 


1.75 


0.629 


0.378 


2.00 


0.782 


0.083 



Table 1: Results of the linear fit ^ logT/3 = W/j^ + n/j of the data in 
Fig. 3. For j3 — > 00, mp — > 1, — > 00. 
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Fig. 1: The critical lengths Ai and A2 as functions of the magnetic field h. 

The solid and dashed lines represent the theoretical prediction (2.3) in the 
limit /9 — > 00; the black circles and the black squares are, respectively, the 
numerical estimates of Ai at = 10 and A2 at /5 = 6. 
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Fig. 2: The solid line represents the average magnetization of the whole 

lattice mo; the other lines represent mi, 1712 iTis and m4. The top figure 
corresponds to the case (3 — 3, h — 0.24 (Ai = 5) and M — 16; the bottom 
one top = 2,h = 0.14 (Ai = 8) and M = 32. 
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23456789 10 

1/h 

Fig. 3: Numerical estimates of logrg//? plotted versus 1/h for different 
values of the inverse temperature (3. Empty circles, black squares, black up- 
ward triangles, black downward triangles and black circles refer respectively 
to (3 = 1.00, 1.25, 1.50, 1.75,2.00. All the results are averages over 60 dif- 
ferent histories for = 32. The solid line is the graph of a linear function 
with slope equal to 1. 
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Fig. 4: Example of a rectangular droplet which is not a local minimum. 
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Fig. 5: Mechanisms of growth peculiar to the Ising model with free boundary conditions. 
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Fig. 6: Protocritical droplet V. 
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